ADP-02-93/T531 



Chiral Extrapolation of Lattice Data for Heavy Baryons 



g ■ X.-H. Guo and A.W. Thomas 

O 
(N 

Q I Department of Physics and Mathematical Physics, 

^ ' and Special Research Center for the Subatomic Structure of Matter, 

^ ■ Adelaide University, SA 5005, Australia 

>• I e-mail: xhguoOphysics. adelaide.edu.au, athomasOphysics. adelaide.edu.au, 

a^ 
cn 

o 

(SI '. Abstract 
O 

• The masses of heavy baryons containing a b quark have been calculated 

numerically in lattice QCD with pion masses which are much larger than 



Oh! 



its physical value. In the present work we extrapolate these lattice data to 
the physical mass of the pion by applying the effective chiral Lagrangian 
^ ' for heavy baryons, which is invariant under chiral symmetry when the light 



quark masses go to zero and heavy quark symmetry when the heavy quark 
masses go to infinity. A phenomenological functional form with three pa- 
rameters, which has the correct behavior in the chiral limit and appropriate 
behavior when the pion mass is large, is proposed to extrapolate the lattice 
data. It is found that the extrapolation deviates noticably from the naive 
linear extrapolation when the pion mass is smaller than about SOOMeV. 
The mass differences between and and between and are also 
presented. Uncertainties arising from both lattice data and our model pa- 
rameters are discussed in detail. We also give a comparision of the results 
in our model with those obtained in the naive linear extrapolations. 



PACS Numbers: 12.39.Fe, 12.39.Hg, 12.38. Gc, 12.40.Yx 



I. Introduction 



The spectrum of some hadrons has been calculated numerically in lattice QCD 
over the past few years. These hadrons include light mesons and baryons [1] , heavy 
mesons [2, 3], and heavy baryons [2, 4]. Using non-relativistic QCD (NRQCD) 
on the lattice [5] for heavy quarks and the tadpole-improved clover action for light 
quarks, the authors of Refs.[2, 3] studied extensively the spectra of heavy mesons and 
heavy baryons (including doubly heavy baryons) . These lattice data were obtained 
in the region where the mass of the pion is much larger than the physical mass of 
the pion. Hence one needs to extrapolate these data to the physical pion mass in 
order to obtain the heavy hadron masses in the real world. Naively, this is done 
by linear extrapolations which are inconsistent with the model independent, non- 
analytic behavior of hadron properties in the chiral limit. In order to overcome 
this problem, pion-hadron loops are included in the study of light hadron properties 
[6, 7, 8, 9]. This yields the correct leading and next-to-leading non-analytic terms 
in the light quark masses and leads to rapid variation at small pion masses. In 
general, lattice data extrapolated to the physical pion mass this way yield quite 
different results from hnear extrapolations. Based on this, we considered previously 
the chiral extrapolation of the lattice data for heavy D and B mesons and discussed 
the important hyperfine splittings [10] . Here we generalize our approach to the case 
of heavy baryons and extrapolate the lattice data for heavy 6-baryons obtained in 
Ref.[2]. 

We work in two opposite limits of quark masses. One is the zero quark mass limit 
while the other is the infinite quark mass limit. When the masses of the light quarks, 
M, d, and s, go to zero the QCD Lagrangian has a chiral SU (3) l x SU (3)^ symmetry, 
which is spontaneously broken into SU{?>)v plus eight Goldstone bosons. When the 
masses of the heavy quarks c and h go to infinity, we have an effective theory, heavy 
quark effective theory (HQET), which is invariant under heavy quark flavor and 
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heavy quark spin transformations, SU{2) / x SU (2)s. Thus the interactions of heavy 
baryons with the hght pseudoscalar mesons should be described by an effective chiral 
Lagrangian for heavy baryons which is invariant under both SU{3)l x SU{3)r and 
SU{2)f X SU{2)s transformations. This chiral Lagrangian will be applied in the 
small pion mass region while we extrapolate the lattice data to the physical pion 
mass. 

The remainder of this paper is organized as follows. In Section II we give a brief 
review of the chiral Lagangian for heavy baryons and including the propagators 
of heavy baryons. In Section III we apply this Lagrangian to calculate pion loop 
contributions to the self-energy of heavy baryons. Then we propose a phenomeno- 
logical functional form with three parameters for extrapolating the lattice data to 
the physical region. In Section IV we use this form to fit the lattice data and give 
numerical results. Finally, Section V contains a summary and discussion. 

II. Chiral perturbation theory for heavy baryons 

When the light quark mass, m^, approaches zero, the QCD Lagrangian possesses 
an SU (3) l x SU (3)^ chiral symmetry. The hght pseudo-Goldstone bosons associated 
with spontaneous breaking of chiral symmetry can be incorporated into a 3x3 matrix 

E = exp(— j. (1) 

where /^r is the pion decay constant, = 132MeV, and M is a matrix which includes 
the eight Goldstones 



M= TT- 

Under SU{3)l x SU{3)r transformations, S is required to transform linearly. 



(2) 



S ^ (3) 
where L e SU{3)l and R e SU{3)r. 



While discussing the interactions of Goldstone bosons with other matter fields 
it is convenient to introduce 

e = Ve. (4) 
Under SU{3)l x SU{3)r transformations 

e ^ m+ = UiR+, (5) 

where the unitary matrix [/ is a complicated nonlinear function of L, R, and the 
Goldstone fields, and is invariant under the parity transformation. 

A heavy baryon is composed of a heavy quark Q {Q — b, or c) and two light 
quarks g^^fi (a. (b) equals 1, 2, 3 for u, d, s quarks, respectively). When the heavy 
quark mass, mg, is much larger than the QCD scale, Aqcd, the light degrees of 
freedom in a heavy baryon become blind to the flavor and spin quantum numbers of 
the heavy quark because of the SU{2) f x SU{2)s symmetries. Therefore, the Hght 
degrees of freedom have good quantum numbers which can be used to classify heavy 
baryons. The angular momentum and parity of the two light quarks may be 
O"*" or 1+, which correspond to SU{3)l+r antitriplet and sextet, respectively. The 
lowest- lying heavy baryons in the 3 representation have spin 1/2, and are denoted 
by fields which destroy these baryons, (T^ — Aq, Ti^2 — Sq). The lowest- lying 
heavy baryons in the 6 representation have spin 1/2 or 3/2, and are denoted by field 
operators S*"* and 5'^*, respectively, where 5'(*)^i'^^'22 _ = '^q\ ^^'^ 

^ ^(*) rp^ transforms under SU{3)l x SU(3)r as 

Ta ^ (6) 
and under heavy quark spin symmetry 

Ta ^ STa, (7) 

where S G SU{2)s. also satisfies 

Ta = ^Ta, (8) 
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where v is the velocity of the heavy baryon. 

It is convenient to combine 5*"^ and 5**"^ into the field 5'^'' [11] 

S;'^^{7, + v,HS^' + S*;\ (9) 

Then under SU{3)l x SU{3)r 

St - U^U'.Sf, (10) 
while under heavy quark spin symmetry 

Sf^SSf. (11) 



S"^^ satisfies the constraints 



Sf = fSf, v^Sf^Q. (12) 



It is convenient to introduce a vector field V^^, 



y^, = lii-'d'i + id>^na,. (13) 
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and an axial- vector field 



K,-\{ed''i-id>^eu. (14) 

Under 5C/(3)l x 5C/(3)r, V^" UV^U+ + Ud>'U+, and A'^ ^ C/A'^C/+. Defining 
the covariant derivative 

(OT)„ = d^^T, - mV't, (15) 

and 

{D^'S^Y^ = d^'Sf + {V^'YSl^ + (l^'')^^r, (16) 
we can show that under 5^/(3)^ x 5t/(3)R, L"^T„ ^ L"^T,,([/+)^ and Di'Sf 

In the limit where light quarks have zero mass and heavy quarks have infinite 
mass, the Lagragian for the strong interactions of heavy baryons with Goldstone 
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pseudoscalar bosons should be invariant under both chiral symmetry and heavy 
quark symmetry. It should also be invariant under Lorentz and parity transforma- 
tions as required in general. The most general form for the Lagragian satisfying 
these requirements is [12] 

^g^{e,,S\A^)\S'^^ + e^'^S^^lAXTal (17) 

where g\ and are coupling constants describing the interactions between heavy 
baryons and Goldstone bosons and AM is the mass difference between sextet and 
antitriplet heavy baryons in the heavy quark limit. As a consequence of heavy quark 
symmetry, g\ and gi are universal for different heavy baryons. Since they contain 
information about the interactions at the quark and gluon level, they cannot be 
fixed from chiral perturbation theory, but should be determined by experiments. 
In the limit mq — > oo, the propagator for Aq is 

i 1 + ^ 
V • p 2 ' 

where p is the residual momentum of the heavy baryon. It can also be shown that 
for Eq, the propagator is 

i 1 + ^ 
vp- AM 2 ' 

and the propagator for Eq is 

1 + ^ -ijgi^u - lln'-yu - lv^^Vu) 1 + ^ 
2 vp-AM 2 

In the limit thq — > oo, there is no mass difference between Eq and Eq. 

In HQET, the leading term which is responsible for a mass difference between 
Eg and Eg is the color-magnetic-moment operator, ^hyO'f^uG'^^hy (where is the 
heavy quark field operator in HQET and G^'^ is the gluon field strength tensor). 
This term is singlet under SU{3)l x SU{3)ji and leads to the following correction 
term to C in Eq.(17): 

i^S^^,a,.S^^\ (18) 
rriQ 
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where a is a constant which also contains interaction information at the quark and 
gluon level, and which is the same for Eq and Eq at the tree level because of 
heavy quark symmetry. When QCD loop corrections are included, a depends on 
niQ logarithmically. 

The term (18) enhances the mass of Eq by a/mq and lowers that of Eq by 
2a /mq. Therefore, the propagators for Eg and Eg become 

i l + i) 



vp-AM+^ 2 ' 



and 



2 VP- AM - ^ 2 ' 

respectively. 

Substituting ^ = exp{iM/ fj,) into Eqs.(13, 14) and making a Taylor expansion 
we obtain the following expressions for and ^4^: 



V, = —^[M,d,M] + 0{M'), 



(19) 



(20) 



Substituting Eq.(9) and Eqs.(19, 20) into Eq.(17) we have the following explicit 
form for the interactions of heavy baryons with Goldstone bosons: 



X 



92 

1 

"73 



1 



(21) 



where O(M^) terms are ignored. 
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Chiral symmetry can be broken explicitly by nonzero light quark masses. This 
leads to the following leading order terms in the explicit chiral symmetry breaking 
masses: 

Ai^„^,(em,e + tm^etS^ + A2^„^,5?IV(m,E+ + Em,) 
+X^TiTa{imqi + tm^tTbT' + A4Tr(f "T,)TrKS+ + Em,), (22) 

where Aj {i = 1,2,3,4) are parameters which are also independent of the heavy 
quark mass in the limit uiq — > oo. 

III. Formalism for the extrapolation of lattice data for 
heavy baryon masses 

Prom the chiral Lagrangian for the interactions of heavy baryons with light Gold- 
stone bosons, Eq.(21), we can calculate pion loop contributions to the heavy baryon 
propagators near the chiral limit - i.e., when the pion mass is not far from the chiral 
hmit. This leads to a dependence of the heavy baryon masses on the pion mass. We 
will concentrate on Eq and Eg, but other heavy baryons can be treated in the same 
way. 

Prom Eq.(21) we can see that there are four diagrams for pion loop corrections 
to the propagator of either Eq or Eg, and three diagrams for Aq. These diagrams 
are shown in Pig. 1 for Eq, in Pig. 2 for Eq, and in Pig. 3 for Aq. It can be easily 
seen that Pig. 1(a), Fig. 2(a) and Fig. 3(a) do not contribute (this is because the 
integrand is of the form k'^f{k^) where k is the momentum of the pion in the loop, 
and f{k'^) is a function of k^) and we will not consider them further. 

Pig. 1(b) arises from the EqttEq vertex. In momentum space it can be expressed 

as 

where p is the residual momentum of the heavy baryon Eg. From Eq.(21), Fig. 1(b) 
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takes the following form: 



9l u v' 



vp- AM + ^ 



(27r)4 [v{p-k)-AM+ ^] - ) ' 



(24) 



where again k is the momentum of the pion in the loop, and is the pion mass. 

As discussed in Ref.[10], the integral 

d^A; ki'k'' 



X 



I 



{2Tr)^ [vk-5]{k'' -ml)' 
where 5 is some constant, can be written as 



(25) 



X^"" = Xi{5)g^'' + X2{5)v''v'', (26) 

where Xi and X2 are Lorentz scalars, which are functions of 5. Obviously, only the 
Xi term contributes in Eq.(24). In the evaluation of Xi, the integration over ko was 
made first by choosing the appropriate contour. Then a cutoff A, which characterizes 
the finite size of the source of the pion, was introduced in the three dimensional inte- 
gration since pion loop contributions are suppressed when the Compton wavelength 
of the pion is smaller than the source of the pion. Since the leading non- analytic 
contribution of these loops is associated with the infrared behavior of the integral, 
it does not depend on the details of the cutoff. In this way, Xi{6) has the following 
expression [10]: 

A + ^/A'^ + mi - 6 m.^ - 5 



X,(5) 



727r2 



12(m^ - 



arctg- 



m1 — S'^ 



arctg 



\ _l_ /a^ + 

+3S{2S^ - 3ml)ln + 36AJA^ + ml + 6{S^ - ml) A + 2A^ 



(27) 



when m^l > 5^; 



X,{6) 



727r2 



6{6^ - ml)'/'\n 



'A + y^A2 + ml -d- _ ^2 



A + ^A2 + ml -5 + ^JS^ - ml 



mt 



m^ 



mt 



+35(2^2 - 3m^)ln ^ + 35A^A2 + ml + 6(5^ - ml) A + 2An , (28) 



m^ 



when <5- In the case where 5 = 0, 



Xi = ^ (^Sm^arctg^ - 3m^A + A' ) . (29) 



Prom Eqs.(23) and (24) we have 



Si = .||xi(AO, (30) 



where Ai = v • p - AM + 

Figs. 1(c), (d) have the same expression as in Eq.(23), except for Si being 
replaced by E2 and E3, respectively. In the same way, we have 

E2=i;^Xi(A2), (31) 

where A2 = f • p - AM —. 

^ f rriQ 

Fig. 1(d) arises from the EqttAq vertex. In this paper we only consider E^*^ 
since lattice data are available for them. For E^, tt^ appears in the pion loop, then 
we have 

^x%X,{A,), (32) 
where /S.^ — v-p. For E°, 7r° appears in the pion loop, and we have 

S? = i^Xi(A3). (33) 

Defining E as the sum of Ei, E2, and E3, the propagator of Eb becomes 

^ 111 iu) 

vp-AM+^-E 2 ' ^ ^ 



where 



= ^ll^i(Ai) + ^^,X^{A,) + i%X,{A,), (35) 



and 



E^° = ^M^l(^l) + ^|L^i(^2) + ^^5^1(^3). (36) 
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Pion loop contributions to the propagator of can be calculated in the same 
way. Fig. 2(b), (c), (d) can be expressed as 

^n^-- irj7 a TT- 9nv - i:l>riv - i:V^v, — — , (37) 



vp-AM-^ 'vp-AM-^ 2 V^'^'^ 3"^"' 3 " V 2 
where i = 1,2,3 for Fig. 2(b), (c), and (d), respectively. After some tedious 
derivations, we obtain 

ni = iMx,(A2), (38) 
n2 = -z-^Xi(Ai), (39) 



and 



n3'^" = 4^i(A3), (40) 

J n 

n3" = i^^i(A3). (41) 



Define 11 as the sum of Hi, 112, and Ha, the propagator of becomes 
1 + ^ -iig^iu - l^f^ju - lv^,Vu) 1 + |^ 



2 vp- AM - ^ - n 



(42) 



where 



and 



n^" = ^M^i(A2) - ^^^i(Ai) + 4^i(A3), (43) 
n"^° = #^i(A2) - ^|,X,(AO + .^X,(A3). (44) 

"/tt "/tt ^/tt 

In the same way, we can calculate pion loop contributions to the propagator of 
A;,. Fig. 3(b) and (c) can be expressed as 

^(-ii^,)^, (45) 
V ■ p V ■ p 

where i = 1, 2 for Fig. 3(b) and (c), respectively. We obtain 

3r72 

Xi = i-|Xi(AO, (46) 



and 



K2^Mx,{A2). (47) 

In 
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If we define K — Ki + K2, then the propagator of A;, becomes 

i 1 + ^ 
vp-K 2 ' 

where 



(48) 



K = t^X,{A,)+z^X^{A,). (49) 

J IT J IT 



After the correction from E is added, the propagator of E;, is proportional to 

1 



(50) 



V ■ p — mo — E(v • p) ' 

where mo is the mass term without E correction. 
The physical mass of E;,, m, is defined by 

[v ■ p - mo - E{v p)]\^.p=m = Q- (51) 

Therefore, to order 0(5'i,5'|) we have 

m = mo + E(t' • p = mo). (52) 

For E^ and A^, E is replaced by 11 and K, respectively, in Eqs. (50-52). For E;,, 
E^, and A^, mo is AM — AM + and 0, respectively. Consequently, the pion 
loop contribution to the mass of E;,, ctsj, has the following expression: 



and 



..2^fA,(0).4x,-^).4A,(AM-|i). (54) 



In the same way. 



and 



» 6/2 6/2 mb fl mb 



= ^M^i(O) - ^^^i(-) + + -)• (56) 

' 6/2 6/2 mfe 2/2 mft 
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For Afe, we have 

^A, = Mx,{-AM + ^)+ i^X,{-/^M - (57) 

In Eqs. (53-57), is given by Eqs. (27-29). 

In order to extrapolate the lattice data from large m^r to the physical value of 
the pion mass, we follow the arguments proposed in Ref.[10] where we dealt with 
heavy mesons. These arguments can be generalized to the case of heavy baryons 
straightforwardly. Eqs. (53-57) arc valid when nij^ is not far away from the chiral limit 
- i.e., when 171^^ < A. As pointed out in Refs.[6, 7, 8, 9, 10], pion loop contributions 
vanish in the limit m,r 00, and the heavy baryon mass becomes proportional to 

when 171.^ becomes large. This behaviour is consistent with lattice simulations. 
Following Refs.[6, 7, 8, 9, 10], we propose the following phenomenological, functional 
form for the extrapolation of lattice data for heavy baryons: 

niB ^ aB + bBml + aB, (58) 

for B = Eb, El or A;,. 

The advantage of fitting the lattice data in this way is that we can guarantee 
that our formalism has both the correct chiral limit behavior and the appropriate 
behavior when uIt^ is large, with only three parameters (a, b, and A) to be determined 
in the fit. 

Chiral symmetry is explicitly broken by the terms in Eq.(22). Substituting 
Eqs.(l, 4, 9) into Eq.(22) we have the following explicit expression: 

2Ai f: [m,a{-SabS^' + s:^s*;')] + 2A2 E ^ {-SabS^' + si^s*;') 

a,b=l a=l a,b=l 

3 3 3 _ 

+2A3 mgafaT'^ + 2A4 rriga ^ f^T", (59) 

a=l a=l a=l 

where we have made a Taylor expansion for ^ and omitted 0(1//^) terms. It can 
be seen that Eq.(59) does not contribute to the mass difference between Eg and 
Eq to order niq. Corrections to this statement are of order mqO{l/ fl), with extra 
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suppression from ruq with respect to the pion loop effects. They will therefore be 
ignored. Eq.(59) may contribute to the mass different between Eg'* and Aq. Such 
effects will be considered to be effectively included in the parameter AM in Eq.(17). 

IV. Extrapolation of lattice data for heavy baryon masses 

The masses of E;,, S^, and were calculated with the aid of NRQCD in quenched 
approximation in Ref.[2]. Since the mass of the heavy quark is much larger than 
A-QCD, it becomes an irrelevant scale for the dynamics inside a heavy hadron and 
is removed from NRQCD. This makes it possible to simulate heavy baryons when 
the lattice spacing is larger than the Compton wavelength of the heavy quark. The 
lattice space used is 1/a — 1.92 GeV. For light quarks the tadpole- improved clover 
action was used which has discretization errors of order a^a. The value of /3 which 
is related to the bare gauge coupling is 6.0 and the lattice size is 16^ x 48. In the 
simulations, three different values for the hopping parameter k, 0.1369, 0.1375, and 
0.13808, were used. The light quark mass is related to k through the definition rUq — 
— 1/kc), with K,c = 0.13917. These three hopping parameters correspond to 
three values of m^: 0.6598GeV2, 0.4833GeV^ and 0.3141GeV2, repectively. 

The heavy baryon masses were calculated for five different values of aM° (M° is 
the bare heavy quark mass): 1.6, 2.0, 2.7, 4.0, 7.0, and 10.0, where the data for the 
last two values are less rehable because of large discretization errors [2] . The best 
estimate for aM°, 2.31, was obtained by matching the lattice data to the mass of the 
B meson. Consequently, in our fit we first extrapolate the lattice data for aM°=1.6, 
2.0, 2.7, 4.0 to aM^ — 2.31. This can be done by hnear extrapolation with respect 
to 1/M° with the form c+ where c and d are constants. This is because aEgim: 
which is the simulation mass in NRQCD and which is related to the heavy baryon 
mass, depends on 1/M° linearly (note that in the case of 6-baryons, 0((1/M°)^) can 
be safely ignored). Then from the data in Table XV of Ref.[2], we obtain the values 
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Table 1: Extrapolated values of aEsim at aM^ — 2.31. 



K 




O-E'sim(Sb) 




0.13690 


0.816(33) 


0.877(28) 


0.889(27) 


0.13750 


0.779(43) 


0.845(32) 


0.856(33) 


0.13808 


0.733(63) 


0.818(40) 


0.827(37) 



of aEsiia for the three hopping parameters at aM° = 2.31, which are shown in Table 
1. In the following, we will extrapolate these values to the physical pion mass with 
the formulas in Eq.(58). 

In our fit we have to determine three parameters in our formalism, (os^^, and 
A in Eq.(58), for example). These parameters are related to AM, ct, gi, and g2, 
which represent interactions at the quark and gluon level and cannot be determined 
from the chiral Lagragian for heavy baryons. In our fit, we treat them as effective 
parameters and assume that their possible slight uiq dependence, which results from 
QCD corrections and I/tuq corrections, has been taken into account effectively in 
this way. 

AM is the mass difference between sextet and antitriplet heavy baryons. Since 
we do not have experimental data for the masses of S^* , we use the data for E^*) 
to determine AM [13]. The spin-averaged mass of E^,*) is |(2ms^ + 4ms*), which 
is bigger than uia^ by 0.213GeV. In our fit, we let AM vary between 0.17GeV and 
0.23GeV, which are given by — rriA^ and ms* — ttia^, respectively. The mass 
difference ms* —m-E^, which is equal to ^ to order l/rric, leads to o; = 0.032GeV^ if 
we choose rric — 0.15GeV. To see the depedence of our fit on a, we let it vary from 
0.025GeV2 to 0.035GeV2. 

The coupling constant g2 can be determined from the decay width for E* — > AcTT, 
which has the following explicit form: 



2 



p 92 
i SS-»Ac7r — 



(m|. -mlf- 2m^(m|. + mj J + 
4mL 



(ms* +mAj^ -ml 



ml. 



(60) 
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Prom ^j:i++^A^T,+ = 18±5GeV, we have g'l = 0.559 ±0.155, while from rs*o^Ac7r- = 
13 ± 5GeV, we have gff = 0.404 ± 0.155. Hence, in our fit we choose the range 
0.249 <gi< 0.714. 

Since E* cannot decay to EcTT, we cannot fix gi from decays. However, gi can be 
related to the matrix of the axial- vector current between sextet heavy baryon states 
where a u ^ d transition is involved. By assuming gj^ — 0.75, which corresponds 
to (y(™icieon _ 25 in neutron /3-decays and using spin-flavor wave functions for heavy 
baryons, the authors in Ref.[14] found that gi = 0.38. Based on this, we let gf vary 
from 0.1 to 0.2 in our fit. 

As discussed in Section HI, the parameter A characterizes the size of the source 
of the pion. In principle, the value of A can be determined by fitting the lattice data. 
However, since A is mainly related to the data at small pion masses and the current 
lattice data are only available at large pion masses, the error in the determination 
of A is very large. The size difference between and E^ is caused by effects of 
order l/m^ which are small. The size difference between E;, and A^ is caused by the 
difference between 0+ and 1+ fight degrees of freedom, which is also the main reason 
for a size difference between N and A. It has been pointed out that the values of A 
for N and A are very close to each other [6]. Hence we expect that the difference 
between the values of A for E^ and A^ should also be small. Since the integrand 
in Xi becomes small near the cutoff A, a small variation in A will only lead to an 
even smaller change in Xi. Based on these arguments, we will ignore the differences 
among the values of A for E^, E^, and A^. To see the dependence of our analysis on 
A, we let A vary between 0.4GeV and 0.6GeV. 

Using the three masses for E^, E^, and A^ in Table 1, we fix the other two 
parameters ( a^^ and bs^ for E^,, for example) in Eq.(58) with the least squares 
fitting method. The values for these two parameters in the cases of E^ and E^^ are 
shown in Table 2, where we choose A = 0.5GeV, a = 0.032GeV2, AM = 0.213GeV, 
gf — 0.15, and 5i| = 0.48. The extrapolated masses for E^, E^, and A^ at the physical 
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Table 2: Fitted parameters, extrapolated masses of , , and and mass 
differences at m^^^^. Numbers in brackets are errors caused by the errors of lattice 
data. 







^b 




a(GeV) 


1.472(0.143) 


1.485(0.187) 


1.290(0.208) 


6(GeV-') 


0.324(0.265) 


0.340(0.326) 


0.439(0.366) 


a(GeV) 




0.0188(0.0026) 




6(GeV-^) 




-0.00166(0.00470) 




m(GeV) 


1.4575(0.1384) 


1.4702(0.1803) 


1.2502(0.2008) 


mjn.± —mj.±{GeV) 




0.0127(0.2272) 




m^T - mA(GeV) 




0.2158(0.2385) 




(ms*± -ms±)*(GeV) 

b h 




0.0180(0.0025) 





pion mass are also shown in this table. The spin-averaged mass is defined as 
i(2ms, + 4ms*). 

With the parameters in Table 2 we obtain the masses of E^, E^, and A^ as a 
function of the pion mass. These results are shown in Figs. 4, 5 and 6, respectively, 
for A = 0.4 and 0.6. The difference between m"^^ and mAj, is plotted in Fig. 7. The 
results of linear extrapolation are also shown in these figures. 

It can be seen from Table 2 that the extrapolated mass difference between E^ 
and E^ has a very large error. This is caused by the uncertainty in the lattice data. 
A better way to obtain the mass difference between E;, and E^ is to extrapolate the 
lattice data for this mass difference which were obtained from ratio fits, since these 
data have much smaller errors. The mass difference between E^ and E^, was 
also given in Ref.[2] for five different values of aM°. We use the data at aM^— 1.6, 
2.0, 2.7, and 4.0 to obtain the value of i^E at aM°=2.31 with the formula 

«AB=;^, (61) 

where e is a constant. Eq.(61) arises from the fact that the mass splitting between 
Efe and E^ is caused by l/mq effects. With the least squares fitting method we 
obtain results for /S.E at aM°=2.31, for different values of k. These are shown in 
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Table 3: Extrapolated values of A£^, the mass difference between and E^, at 
aM^ = 2.31 using Eq.(61). 



K 


0.13690 


0.13750 


0.13808 


aAE 


0.0093(5) 


0.0095(7) 


0.0096(7) 


AE 


0.0179(10) 


0.0182(13) 


0.0184(14) 



Table 3. 

In order to extrapolate the values in Table 3 to the physical mass of the pion, 
we use the following formula: 

ms* - = a + bml + cts* - ctsj. (62) 

With A = 0.5GeV, a = 0.032GeV2, AM = 0.213GeV, gl = 0.15, and = 0.48, 
we obtain a — 0.0188(26), b — —0.00166(470), and the extrapolated mass diference 
between E^ and E^^ m-^*± — m-^± = 0.0180(25), which is hsted in Table 2 as 

b b 

{'mj^*± — mj^±)*. In Fig. 8, we show mj..± — mj;± obtained in this way as a function 

6 6 6 6 

of the pion mass. For comparison, in Fig. 9 we plot the result for 'm-^*± — "^^i which 

6 6 

is obtained from Eq.(58), although it has very large errors. Because these errors are 
so large, the extrapolation from Eq. (58) is consistent with the extrapolation directly 
from the lattice data for the mass difference between E^ and E^. 

In addition to the uncertainties which are caused by the errors in the lattice 
data, the fitted results can also vary a little in the range of the parameters a, AM, 
g\, g2, and A. In Table 4 we list these uncertainties. 

In the naive linear extrapolations pion loop corrections are ignored. Hence the 
results do not depend on the parameters a, AM, g\, gf|, and A. In Table 5 we 
list the results of linear extrapolations for comparison. We note that there is no 
difference between the results for E^*^^ and E^*^*^ in the hnear extrapolations. 

Comparing the uncertainties listed in Table 2 and Table 4 we can see clearly that 
the main uncertainties in our fit are caused by the errors in the lattice data. In fact. 
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Table 4: Uncertainties for the extrapolated quantities for and which are 
caused by the uncertainties associated with parameters in the fitting function. 



Quantities 


a 


AM 




9l 


A 


niyo 


0.02% 


0.3% 


0.06% 


0.8% 


0.9% 




0.007% 


0.3% 


0.04% 


0.9% 


1.0% 




0.0% 


0.2% 


0.0% 


1.7% 


1.8% 


ms*o — m^o 

b h 


3.1% 


0.8% 


1.6% 


12.6% 


10.2% 




0.0% 


3.4% 


0.3% 


3.7% 


9.3% 


(ms*o - m^o)* 


2.2% 


0.0% 


1.1% 


7.7% 


1.7% 



Table 5: Fitted parameters, extrapolated masses of E^, E^, and A^ and mass differ- 
ences at mP^^^ for hnear extrapolations. Numbers in brackets are the errors caused 

by the errors of the lattice data. 







E* 


Aft 


a(GeV) 


1.465(0.143) 


1.479(0.187) 


1.63(0.208) 


6(GeV^^) 


0.330(0.265) 


0.346(0.326) 


0.460(0.366) 


a(GeV) 




0.0190(0.0025) 




6(GeV-^) 




-0.00172(0.00470) 




m(GeV) 


1.4714(0.1384) 


1.4854(0.1803) 


1.2724(0.2008) 


/;/x,-.± — /;k-.± (GeV) 




0.0110(0.2272) 




ni^^ - nixiGcV) 




0.2081(0.2385) 




(ms*± -ms±)*(GeV) 

b h 




0.0190(0.0025) 
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the errors of lattice data for heavy baryons are much larger than those for heavy 
mesons [3]. Indeed, the uncertainties in the extrapolated heavy baryon masses are 
about one order larger than those in the case of heavy mesons. However, because 
of the small errors in the lattice data for the mass splitting between and E^^ 
the extrapolated mass difference at the physical pion mass also has a smaller error, 
about 28%. 

Prom Figures 4 to 9 we see that when the pion mass is smaller than about SOOMeV 
the extrapolations begin to deviate from linear behavior. This is because the pion 
loop corrections begin to affect the extrapolations around this point. As the pion 
mass becomes smaller and smaller, pion loop corrections become more and more 
important. For the masses of E^, E^^, A^, and the mass difference between E^ and 
E^^ the extrapolated values are smaller than those obtained by linear extrapolation. 
For the difference between the spin-averaged mass of eJ*^^ and the mass of Afe, the 
extrapolated value is larger than that obtained by linear extrapolation. We have 
checked that such behaviour is independent of the uncertainties in the parameters 
in our model. 

Comparing the results in the naive linear extrapolations with those with pion 
loop corrections being included we find that the difference between them is much 
smaller than that found in the case of mesons. For example, the splitting between 
E^ and E^^ is about 5% smaller if pion loop effects are taken into account, while 
the hyperfine splitting in the case of B mesons is about 20% smaller when pion 
loop effects are taken into account [10]. For the masses of E^, E^^, and A^, the 
extrapolated values with pion loop effects being included are about only 1% smaller 
than those in linear extrapolations, while for B and B* the corresponding number 
is about 3% [10]. Hence in the case of heavy baryons, the linear extrapolation is a 
better approximation than in the case of heavy mesons. 

For E° and E^° we should use Eqs.(54, 56) in the extrapolation of lattice data. 
Repeating the same procedure as that for E^ and E^^ we find that, in addition 
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Table 6: Fitted parameters, extrapolated masses of and and mass difference 
at mP^^^. Numbers in brackets are errors caused by the errors in the lattice data. 





yO 


^b 


a(GeV) 


1.469(0.143) 


1.482(0.187) 


6(GeV-^) 


0.327(0.265) 


0.342(0.326) 


a(GeV) 


0.0187(0.0026) 




6(GeV-^) 


-0.00151(0.00470) 




m(GeV) 


1.4638(0.1384) 


1.4774(0.1803) 


m^'i — iTij.± (GeV) 


0.0135(0.2272) 




m^^l -mA(GeV) 


0.2226(0.2385) 




(ms*± -ms±)*(GeV) 

b b 


0.0187(0.0025) 





Table 7: Uncertainties for the extrapolated quantities for E" and E^ which are 
caused by the uncertainties in the parameters of the fitting function. 



Quantities 


a 


AM 


9! 


9l 


A 




o.orx 


0.2% 


0.06% 


0.1% 


0.1% 




0.007% 


0.1% 


0.04% 


0.5% 


0.4% 




0.0% 


0.2% 


0.0% 


1.7% 


1.8% 


m^*± — mv± 


1.5% 


0.0% 


2.2% 


5.9% 


1.5% 




0.0% 


2.3% 


0.3% 


6.6% 


7.9% 


(ms*± - m^±)* 

b b 


1.6% 


0.0% 


1.1% 


0.7% 


0.5% 



to some minor changes in numerical results, all the quantatitive results remain the 
same. In Tables 6 and 7 we hst our numerical results for E° and E^°. Comparing 
the results in Table 6 with those in Tables 4 and 5 we can see that the naive hnear 
extrapolations work even better in the case of E^ and E^° than in the case of E^ 
and E^±. 

V. Summary and discussion 

The masses of heavy baryons E^, E^, A;,, and the mass difference between E;, and 
E^ have been calculated numerically in lattice QCD with unphysical pion masses 
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which are larger than about 560MeV. In order to extrapolate these data to the 
physical mass of the pion in a consistent way, we included pion loop effects on the 
heavy baryon masses by applying the effective chiral Lagrangian for heavy baryons 
when the pion mass is smaller than the inverse radii of heavy baryons. This chiral 
Lagrangian is invariant under both chiral symmetry (when the light quark masses go 
to zero) and heavy quark symmetry (when the heavy quark masses go to infinity). 
In order to study mass difference between and E^, we took the color- magnetic- 
moment operator at order l/mq in HQET into account since this operator is the 
leading one to cause splitting between E^ and E^. When ttItt becomes large, lattice 
data show that heavy baryon masses depend on m\ linearly in the range of interest. 
Based on these considerations, we proposed a phenomenological functional form to 
extrapolate the lattice data. 

The advantage of our formalism is that it has the correct chiral limit behavior 
as well as the appropriate behavior when is large and that there are only three 
parameters to be determined in the fit to lattice data. It is found that when the 
pion mass is smaller than about SOOMeV the extrapolations begin to deviate from 
the naive linear extrapolations. However, the differences between the extrapolations 
with and without pion loop effects being included is smaller than those in the case of 
heavy mesons. Hence the linear extrapolation is a better approximation in the case 
of heavy baryons. We carefully analysed uncertainties in our extrapolations which 
are caused by both lattice data errors and uncertainties in several parameters in our 
model and found that the main uncertainties are caused by the errors of the lattice 
data. By directly extrapolating the lattice data for ms* — ^n^i which has much 
smaller errors, we found that the extrapolated mass difference between E^ and E^^ 
at the physical mass of the pion is 18.0 MeV, with an uncertainty of 28% caused by 
lattice data errors. For E^ and E^" this difference is 18.7MeV with 26% uncertainty 
from lattice data errors. The uncertainties associated with the parameters in our 
model are at most a few percent. For the difference between the spin-averaged mass 
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of and the mass of Aft, the extrapolated value has a very large error. This 

needs to be improved when the lattice data become more accurate. Furthermore, 
we should bear in mind that our extrapolations are based on the lattice data in 
the quenched approximation. Prom our experience in the cases of light and heavy 
mesons [15, 10], the quenched approximation may affect the mass splitting between 
Eft and E^. In addition, the lattice results for m^* — msj, may be sensitive to both 
the coefficient of the a ■ B term in NRQCD [3] and the clover coefficient in the clover 
action for light quarks. This may also influence the lattice data and consequently 
affect our extrapolations. 
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Figure Captions 



Fig. 1 Pion loop corrections to the propagator of SU (3) sextet heavy baryons with 
spin-1/2, where Sq^ represent spin-l/2(3/2) SU{3) sextet heavy baryons with heavy 
quark Q and Tq represents SU (3) antitriplet heavy baryons. 

Fig. 2 Pion loop corrections to the propagator of SU (3) sextet heavy baryons with 
spin-3/2. Same notation as in Fig. 1. 

Fig. 3 Pion loop corrections to the propagator of SU{3) antitriplet heavy baryons. 
Same notation as in Fig. 1. 

Fig. 4 Phenomenological fits to the lattice data for the masses of as a function 
of the pion mass. The solid (dashed) line corresponds to A = 0.4(0. 6)GeV and the 
dotted line represents the fit using a linear extrapolation. 

Fig. 5 Phenomenological fits to the lattice data for the masses of as a function 
of the pion mass. Same notation as in Fig. 4. 

Fig. 6 Phenomenological fits to the lattice data for the masses of A;, as a function 
of the pion mass. Same notation as in Fig. 4. 

Fig. 7 Difference between the spin-averaged mass of and the mass of A5 as a 

function of the pion mass, which is obtained from Figs. 4, 5 and 6 (same notation as 
in Fig. 4). 

Fig. 8 Phenomenological fits to the lattice data for the mass difference between 
and E^^ as a function of the pion mass (same notation as in Fig. 4). 

Fig. 9 mv*± — my;± as a function of the pion mass, which is obtained from Figs. 4 

b b 

and 5. The large errors of lattice data are not shown (same notation as in Fig. 4). 
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